On Estimation of Hurst Scaling Exponent 
through Discrete Wavelets 
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Abstract 

We study the scaling behavior of the fluctuations, as extracted through wavelet 
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coefficients based on discrete wavelets. The analysis is carried out on a variety of 

physical data sets, as well as Gaussian white noise and binomial multi-fractal model 
fS) \ time series and the results are compared with continuous wavelet based average 

wavelet coefficient method. It is found that high-pass coefficients of wavelets, be- 
longing to the Daubechies family are quite good in estimating the true power in the 
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fluctuations in a non-stationary time series. Hence, the fluctuation functions based 
on discrete wavelet coefficients find the Hurst scaling exponents accurately. 
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1 Introduction 



Fractals exhibiting self-similar behavior, are ubiquitous in nature (|lj). They 
manifest in areas ranging from financial markets to natural sciences. Several 
techniques have been developed to study the correlations and scaling prop- 
erties of time series exhibiting self-similar behavior; some of these data sets 
are non- stationary in character. Various methods like R/S analysis (|2|), struc- 
ture function method (J3|), wavelet transform modulus maxima (|4j), detrended 
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fluctuation analysis and its variants p; la; 13; |s|; [9|; llO; llll; LUl; UM, LljJ), aver- 



age wavelet coefficient method (Il5l ) and a recently developed discrete wavelet 
based approach by the present authors, have been employed for the charac- 
terization of fluctuations (1 16c [lj 



Wavelets, through their multi-resolution and localization abilities, are well 
suited for extracting fluctuations at various scales from local trends over ap- 



propriate window sizes (Il8l;ll9l). The nature of the fluctuations extracted partly 
depend on the choice of the wavelets, which are designed to have properties 
useful for a desired analysis. For example, the Daubechies family of wavelets 
satisfy vanishing moment conditions, which make them ideal to separate poly- 
nomial trends in a data set. In a number of wavelet based approach for char- 
acterizing self-similar data, wavelet high-pass coefficients are used in finding 
the fluctuation function or other related quantities. For example, in the aver- 
age wavelet coefficient method, continuous wavelet transform is used to find 
the variations through the corresponding wavelet coefficients, from which the 
Hurst scaling exponent is computed. 



The goal of the present article is to analyze the self-similar properties of the 
fluctuations, as extracted through discrete wavelet coefficients, for the pur- 
pose of checking its efficacy vis-a-vis the continuous wavelet based method. 
The study is carried out on Gaussian white noise and binomial multifractal 
model time series, as well as a number of physical data sets. It is found that 
the fluctuation function based on discrete wavelet coefficients is accurate in 
estimating Hurst scaling exponent. It is worth emphasizing that continuous 
wavelets, because of their over completeness tend to overestimate the power. 
In comparison, the discrete wavelets provide a complete orthonormal basis, 
ensuring that fluctuations are independent at each level. 



The paper is organized as follows. In the following section, we study the nature 
of the fluctuations extracted through the wavelet coefficients. In Sec. Ill, we 
then proceed to the detailed analysis of fluctuations of Gaussian white noise 
and binomial multi-fractal model and compare the discrete wavelet based 
result with average wavelet coefficient method. Subsequently, we carry out 
analysis of the fluctuations in the data of observed ionization current and 
potentials in tokamak plasma, time series constructed from random matrix 
ensembles and the financial data sets belonging to NASDAQ and Bombay 
stock exchange (BSE) composite indices. Finally, we conclude after summa- 
rizing our findings and giving future directions of work. 



2 Fluctuations in the wavelet domain 



In discrete wavelet transform it is well-known that, a given signal belonging to 
L 2 space can be represented in a nested vector space spanned by the scaling 
functions alone. This basic requirement of multi-resolution analysis (MRA) 



can be formally written as (j20l ) 



... C V-2 C V-\ C z/_o C1/1C v 2 ... C L 2 , (1) 



with z/.oo = and v^ = L 2 . This provides a successive approximation of 
a given signal in terms of low-pass or approximation coefficients. It is clear 
that, the space that contains high resolution signals will also contain signals 
of lower resolution. The signal or time series can be approximated at a level 
of ones choice, for use in finding the local trend over a desired window. The 
fluctuations can then be obtained by subtracting the above trend from the 
signal. We have followed this approach earlier for extracting the fluctuations, 
by elimination of local polynomial trends through the Daubechies wavelets 

HQ. 

Wavelets also provide a decomposition of a signal in terms of wavelet coeffi- 
cients and one low-pass coefficient: 

L 2 = ... © W^ 2 © W_i © W © Wi © W 2 - (2) 



and 

W^oo © ... © W-! = Uq. (3) 



Wavelet coefficients represent variations of the signal at different scales. For 
example, level one coefficients capture the highest frequency components, cor- 
responding to variations at highest resolution and other wavelet coefficients 
represent variations at progressively higher scales or lower resolutions. As men- 
tioned earlier, these coefficients can differ significantly from the true fluctua- 
tions in the data sets. Below, we explore this aspect through the estimation 
of Hurst exponents. 

Let Xt (t = 1, 2, 3, ..., N) be the time series of length N. First one determines 
the "profile" (say Y(i)), which is cumulative sum of series after subtracting 
the mean. 



n*) = Et** -<*>]> i = l,-;N. (4) 



Next, we obtain the statistics of scale dependence by transforming the profile 
of the time series into wavelet space, the coefficients of wavelets at various 
scales s are used to determine the fluctuation function. The high frequency 
details are captured by the lower scale wavelet coefficients and the higher 
scales capture the low frequency details. By convolving the discrete wavelet 
transform over the given time series, the wavelet coefficients are obtained for 
various scales: 

AT-l 

w j * = v' 2 j2 Y i Hvt-k). (5) 

j=0 



Here 'j' is the scaling index and k represents the translation variable. Since 
discrete wavelet transform satisfies orthogonality condition, it can provide the 
information of time series at various scales unambiguously. Performing wavelet 
transform using Daubechies basis, the polynomial trends in the time series are 
eliminated. In the analysis carried out below we make use of the Daubechies-4 
wavelets. As has been observed earlier, small fluctuations are least affected by 
this basis and straight line trends (akin to a linear fit) are removed through 
the use of this wavelet (Hi 



The wavelet power is calculated by summing the squares of the coefficient 
values for each level: 



m = E w&- (6) 

fc=0 



To characterize the time series, the fluctuation function F(s) at a level s is 
obtained from the cumulative power spectrum: 



F(s) 



E^O') 



1/2 

(7) 



The scaling behavior is then obtained through, 

F(s) ~ s H . 



Here H is the Hurst scaling exponent, which can be obtained from the slope 
of the log-log plot of F(s) vs scales s. It is well known that Hurst exponent is 
one of the fractal measures, which varies from < H < 1. For persistent time 
series H > 0.5 and H = 0.5 uncorrelated series. H < 0.5 for anti-persistent 
time series. 
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Fig. 1. Time series of (a) ion saturation current (IC), (b) floating potential (FP), 
6 mm inside the main plasma and (c) ion saturation current (ISC), when the probe 
is in the limiter shadow. Each time series is of approx. 24, 000 data points. 




Fig. 2. Time series of (a) NASDAQ composite index for a period of 20 years, 
starting from ll-Oct-1984 to 24-Nov-2004, and (b) BSE sensex index, over a period 
of 15 years, starting from 2-Jan-1991 to 12-May-2005. 
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Fig. 3. Wavelet coefficients at various levels for NASDAQ composite data. 
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Fig. 4. Wavelet coefficients at various levels for tokamak plasma data involving ion 
saturation current from top to bottom respectively. 

We have analyzed three sets of experimentally observed time series of variables 
in ohmically heated edge plasma in Aditya tokamak (I2ll ). The time series are 
i) ion saturation current, ii) ion saturation current when the probe is in the 
limiter shadow, and iii) floating potential, 6mm inside the main plasma. Each 
time series has about 24,000 data points sampled at 1MHZ (I22l ). These are 
shown in Fig. 1. The study of fluctuations play an important role in our 
understanding of turbulent transport of particles and heat in the plasma. In 
Fig. 2, we show financial time series of NASDAQ composite index and BSE 
sensex index values. Wavelet coefficients at various scales have been displayed 
in Fig. 3 and Fig. 4. In Figs. 5 and 6, we have shown F(s) versus s for the 
time series of three experimental data sets and financial stock market data, 
respectively. The scaling exponent H, for all the three experimental time series 
as well as financial data sets shows long range correlations (H > 0.5). We have 
also analyzed the discrete time series obtained from random matrix ensembles 
corresponding to Gaussian orthogonal ensemble (GOE), Gaussian symplectic 
ensemble (GSE) and Gaussian unitary ensemble (GUE). These show long 
range anti-correlation behaviors H < 0.5. Gaussian diagonal ensemble (GDE) 
shows uncorrelated behavior, H = 0.5. It is worth mentioning that, we have 
followed the recent approach of Ref. (|23f ) . for converting the random matrix 
ensemble data to discrete time series (|24j). In Table-I and II, Hurst exponents 
of various data sets are given. These results agree with our previous discrete 
wavelet based approach. 



3 Conclusion 



In conclusion, one needs to be careful in using wavelets for estimating the 
Hurst exponent. The continuous wavelets provide an over complete basis, be- 
cause of which the coefficients are not independent, they overestimate the 
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Fig. 5. The log-log plot of fluctuation function F(s) vs s, for the time series of 
NASDAQ composite index and BSE sensex index values for returns and shuffled 
returns. One clearly sees long range correlation behavior. 




Fig. 6. The log-log plot of fluctuation function F(s) vs s, for the time series of 
tokamak plasma data. For larger window sizes one observes long-range correlations. 

power in the fluctuations. In comparison the discrete wavelet is based on com- 
plete orthonormal basis, ensuring that fluctuations are independent at each 
level. Hence this discrete wavelet approach yields correct values for the Hurst 
exponent. 
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Hurst (H) 


NASDAQ - returns 
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IC 
ISC 
FP 
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0.554 
0.549 


GOE 
GDE 
GSE 

GUE 
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0.143 
0.107 



Table 1 

Computed Hurst exponent for various data sets, involving financial, tokamak plasma 

and random matrix energy fluctuations. 
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H a 


H d 


H c 


BMF 

WGN 
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0.5077 


0.8433 
0.5163 



Table 2 

Hurst scaling exponent computed analytically (H a ), through discrete wavelet co- 
efficient method (i^)and continuous average wavelet coefficient method (H c ) for 
both White Gaussian Noise(WGN) and Binomial Multi-Fractal Model (BMF) time 
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